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The electronic and transport properties of graphene modulated by magnetic barrier arrays are 
derived for finite temperature. Prominent conductance gaps, originating from quantum interference 
effects are found in the periodic array case. When a structural defect is inserted in the array, 
sharp defect modes of high conductance appear within the conductance gaps. These modes can 
be shifted by local doping in the defect region resulting into sensing of the chemical molecules 
that adhere on the graphene sheet. In general it is found that sensitivity is strongly dependent on 
temperature due to smoothing out of the defect-induced peaks and transport gaps. This temperature 
dependence, however, offers the added capability for sub-mK temperature sensing resolution, and 
thus an opportunity towards ultra-sensitive combined electrochemical-calorimetric sensing. 


I. INTRODUCTION 

Advances in fabrication techniques since the isola¬ 
tion of monolayer graphene^ have allowed the realiza¬ 
tion of a variety of graphene-based applications^!^. In 
terms of electronic applications, graphene’s high carrier 
mobility^!^ has enabled its application as an electrical 
conducting channel^^— , however, the lack of strong con¬ 
ductance modulation^ due to Klein tunnelingiS^— has 
limited its utilization for other graphene-based electronic 
devices. The modulation of charge carriers in graphene 
is still a significant issue of ongoing research for graphene 
nanoelectronicsi^ii^. 

A conductance modulation in graphene can be in¬ 
troduced by utilizing various mechanisms that open 
an electronic band gap, such as the interaction with 
susbtratesi^^— , elastic strainiSr^, or finite-size effects 
in graphene nanoribbons22i2^. Aside from the band gap 
opening, however, these methods also introduce unavoid¬ 
able disorders, such as charge impurities or structural 
imperfections by substrates^ir— and strong backscatter- 
ing by rough edges of graphene nanoribbons^SrJS, that 
deteriorate the transport properties of graphene. An al¬ 
ternative mechanism for band-gap opening, that in prin¬ 
ciple does not produce any disorder effects in graphene, is 
the use of inhomogeneous magnetic fields. Particularly, 
since inhomogeneous magnetic fields can tune the Dirac 
fermion transport in graphene mimicking a typical po¬ 
tential barrier for Dirac fermion o^^’^^ , one can produce 
magnetic confinement effects for graphene carriers22i2i. 
These effects have motivated research on various struc¬ 
tures, e.g., magnetically defined quantum dots^i^^, Dirac 
fermion waveguide a^'^’^^ , and superlattice a ^ i ^. 

Besides nanoelectronics, another important aspect of 
graphene applications is gas sensin g^^d° . The operation 
principle of a graphene gas sensor relies on measuring the 
modulated transport properties of graphene that are in¬ 
duced by the adhesion of gas molecules on the graphene 
surface. In particular, conductance changes arise from 
the injection of charged carriers into graphene from ad¬ 
hered gas molecules such as O 2 , NO 2 or NH 3 , which play 


roles of electron donors or acceptors^—. It has been 
shown that highly sensitive graphene-based gas sensing 
capable of detecting individual molecule adhesion^ is 
possible, but a large array of sensors might still be re¬ 
quired in order to enlarge the exposed area of graphene 
and thus increase the chances for molecule adhesion at 
shot-time exposures and minute concentrations. Alterna¬ 
tively, here we explore ways of increasing the minimum 
exposed area required for single molecule detection with¬ 
out making arrays, which is significant for minimizing 
device sizes and detection speeds. 

Here, we theoretically study an alternative possibility 
for achieving highly sensitive gas sensing, which is to uti¬ 
lize the adsorbate-induced electrochemical doping effects 
in order to modulate the tunneling resonances (TR) and 
transport band gaps (TBG) emerging in magnetic barrier 
arrays in graphene. We show that in such graphene sen¬ 
sors the presence of adsorbates introduces strong changes 
into graphene’s electrical conductance. We take into ac¬ 
count two structures for the conductance modulation, as 
illustrated by Fig. [TJ First, we consider a periodic array 
of magnetic barriers in which the inter-barrier regions 
are exposed to electrochemical adsorbates. In the second 
case, we consider a periodic magnetic barrier array with 
a structural defect, where electrochemical doping is in¬ 
duced only in the defect region. We investigate the ballis¬ 
tic conductance and its modulation by doping effects for 
both structures, in different temperature ranges. Specif¬ 
ically, in the zero-temperature limit, the presence of a 
structural defect leads to sharp resonance peaks in the 
conductance spectra, that can be sensitively modulated 
by local doping in the defect region. In finite-temperature 
cases, on the other hand, the overall sensitivity is re¬ 
duced because of thermal smoothing, but still remains 
significant within the transport band gaps when a peri¬ 
odic magnetic barrier array with electrochemical doping 
in-between the barriers is considered. By the same token, 
however, the strong temperature dependence which is re¬ 
garded as a weak point in terms of gas sensing, can be 
used for temperature sensing. The highest sensing ability 
is expected at lower temperatures, and yields a promising 
sensing platform for applications in a low-temperature 
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FIG. 1. (a) and (b) Schematic views of magnetic barriers considered in the zero-temperature limit and finite-temperature 

regime. A series of ferromagnetic (FM) stripes (red blocks) are placed upon graphene sheet, and adsorbates (blue spheres) can 
be adhered in specific regions, (c) and (d) Magnetic vector potentials Ay and their corresponding magnetic field Bz profiles 
for the magnetic structures (a) and (b), respectively. Shaded boxes implies adsorbate-induced local doping potentials. The 
magnetic structures are characterized by the barrier width W, the inter-barrier distance D, and the width of the defect region 
Wi. 


environment. 

This manuscript is organized in the following manner: 
in Sec. m we present the model Hamiltonian employed 
to represent multiple magnetic barriers, and we explain 
the transfer matrix formalism which is valid in the bal¬ 
listic regime at finite temperatures. In Sec. uni we dis¬ 
cuss the features of the defect-induced tunneling reso¬ 
nances in the transmission spectra. Next, in Sec. m we 
show the calculated results for the sensing effect by lo¬ 
cal doping in the zero-temperature limit, and in Sec. 0 
we discuss the large conductance modulation beyond the 
zero-temperature limit and its temperature dependence. 
Finally, section IVII contains the conclusions and a sum¬ 
mary of our results. 


II. MODEL HAMILTONIAN AND FORMALISM 


where vp ~ 10® m/s is the Fermi velocity of Dirac 
fermions, and a = {cri,a 2 ) are Pauli matrices acting on 
sublattices of graphene. For simplicity, the magnetic bar¬ 
rier is characterized by a rectangular profile of the vector 
potential: 

A (x) = BIb [0 (a;) - 0 (a; - W)] y, (2) 


where B is the magnetic field strength, Ib = ^/h/eB 
the characteristic magnetic length, 0 (a:) is the Heavi¬ 
side step function and W is the width of the magnetic 
barrier. The corresponding magnetic fields are given by 
B = V X A, which corresponds to delta-function-like 
spikes of opposite sign at the two edges of the magnetic 
barrier. A general expression of the solution in the three 
different regions a;<0, 0 <a:<FF and IF < a; is written 
as^: 


We consider a graphene sheet with a periodic array 
of magnetic barriers. Magnetic barriers have been ex¬ 
perimentally realized by ferromagnetic stripes on top of 
a graphene sheet^. To describe Dirac fermion ballistic 
transport through the system we will use the transfer 
matrix formalism. Starting with the simplest case of a 
single magnetic barrier along the y-direction, the Dirac 
Hamiltonian reads^ 


i; (X) = ) , 

(3) 

where j = 1,2, or 3 represents different regions. The 
solution can be written in the matrix form: 


p + eA (x) 


+ U{x), 


( 1 ) 


(x) 


Qj (x) 



H = ■ 


(4) 
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with the matrices Qi = Q 3 = Q and Q 2 = Q' dehned 
as: 


Q(a:) 

Q' (x) 




^—ikx 


_gg-zfc*Xg-Z(/i 


gik'^x 
g! ^ik'^x 



(5) 


where 


k 


X — 


/ 


Z-U 

hvF 


2 


1.2 

Ky, 



s = sgn(e - U ), 



4>' = tan ^ 


s' = sgn (e). 


^ /Cy + CAyjfl^ 

( 6 ) 


We take into account the low-lying excitation of Dirac 
fermions in graphene, below 100 meV from the charge 
neutral point, in order to remain within the linear-band 
approximation governed by the Dirac equation. The 
magnetic barrier height is found from the definition of 
k'^ at normal incidence to be Eb = vpeAy. To avoid hav¬ 
ing negligible tunneling of Dirac fermions through the 
magnetic barrier, the magnetic barrier width should not 
exceed a few tens nm. Additionally, the total lateral size 
of a multiple magnetic barrier array should not be over 
a few microns so we could still be in ballistic transport 
regime even at room-temperatures^. 

In order to obtain the transmission and reflection prob¬ 
abilities through a single magnetic barrier, we calculate 
the undetermined coefficients at and hi through the ap¬ 
plication of the boundary conditions, i.e. that wavefunc- 
tions should be continuous at the interfaces x = 0 and 
X = W. We note that in this study we do not need to take 
into account the Zeeman energy a ■ gfisB, where g is the 
gyromagnetic factor for Dirac fermions in graphene and 
/rs is the Bohr magneton, because the anti-symmetric 
magnetic field profiles of each single magnetic barrier 
yield no spin-dependent transport phenomenal^. The 
wavefunction continuity provides the following equations: 

Q(0)(j;)=Q'(0)(“j). 


Alternatively, the transfer matrix could also be derived 
using optical analogies, as shown in Appendix!^ We note 
here that since we are interested in ballistic transport of 
Dirac fermions through the structures, no disorder effects 
that could lead to energy dissipation by inelastic scatter¬ 
ing were included in the our formalism. 

The transmission probability T is obtained from the 
transfer matrix by T = \tf. Note that T is not only a 
function of the energy but also of the incident angle of 
the incoming Dirac fermions. In the ballistic regime and 
at finite temperature, the conductance through a two- 
dimensional system is obtained by the weighted average 
of the transmission function over the incident angle, in 
accordance to the Landauer-Biittiker formalism^: 

= / er(e,0)cos<^ 

TT^h^^VF J -00 J-zv/2 

where Ly is the system size in the transverse 
(y) direction, </) is the propagating direction of 
incoming Dirac fermions, and f{e,EF,T) = 
{1 + exp[{e — Ef)/ kBT]}~^ is Fermi-Dirac distri¬ 
bution with given Fermi energy Ef and temperature T. 
In the zero temperature limit the conductance formula 
is simplified to: 


Ap'^T ^+7r/2 

G = 2 J ^ COS 

n^h^VF J-n/2 

( 11 ) 

Next, we consider a monolayer graphene sheet in the 
presence of a series of magnetic barriers. The models 
considered in this study are shown in Fig.[Tl where we as¬ 
sume periodically arranged barriers with an inter-barrier 
distance D. Note that we set an alternating configura¬ 
tion of the magnetic barriers, i.e., alternating signs of 
the vector potentials. The transfer matrix formalism for 
multiple magnetic barriers can be extended from the case 
of the single magnetic barrier: 

N 

T = nT*, ( 12 ) 

2=0 


Q'(K')(“^^)=Q(ir)(;=). (7) 

By combining these, we get an equation connecting the 
coefficients the incoming and outgoing solutions: 



where N is the number of magnetic barriers in use, and 

r ci-Hx, + w)q,'+{xi + w) 

_ I xQ(r^ (Xi) Q (Xi), Hs odd 

* I Q-~^{x^EW)q,'_{x^+W) 

I ^ Q'r' (x,) Q (xi), i is even 

(13) 


where T is the transfer matri x'^^d^ . 


T = Q-i (IF) Q' (IF) Q'-i (0) Q (0), 



where Xi = i (IF + D) is the position of f-th magnetic 
barrier and Q± is for upward or downward magnetic 
barriers. In the following sections we assume 20 mag¬ 
netic barriers periodically arranged with IF = 50 nm, 
D = 50 nm, Eb = VF^Ay = 20 meV and Wd = 100 nm. 
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III. TRANSMISSION PROBABILITY 
THROUGH MULTIPLE MAGNETIC BARRIERS 

In this section, we discuss defect-induced resonances 
and local-doping effects in the transmission probability of 
Dirac fermions through multiple magnetic barriers. The 
existence of defect-induced resonances is best demon¬ 
strated by comparing the transmission spectra through 
a periodic array of magnetic barriers without and with 
a structural defect, as illustrated in Fig. [TJa) and (b), 
respectively, and plotted for normal incidence in Fig. [51 


assume local doping applied to the defect region only, 
as respectively illustrated in Fig. |TJa) and (b). Doping 
is introduced in our model as a rectangular electrostatic 
potential barrier, characterized by the barrier height Ud- 



Energy (meV) 


FIG. 2. Transmission probabilities for the normal incidence 
versns Fermi energy through multiple magnetic barriers with 
and without a structural defect, in linear and logarithmic 
scales. Shaded regions represent transporting gaps. 

It is shown that the periodic arrangement of magnetic 
barriers leads to the existence of transport gaps (TG) in 
specific energy ranges, indicated by shaded regions in Fig. 
121 The emergence of TGs comes from the combination 
of strong back scattering by the magnetic potentials and 
quantum interference effects in the periodic structure, 
similarly to the Kronig-Penny model. (The existence of 
the TG is best understood by the band structure of the 
infinite magnetic barrier array, shown in Appendix|Bl) At 
energies below the magnetic barrier height (e < 20 meV), 
only resonant tunnelings through the multiple magnetic 
barrier are allowed, with the number of resonances be¬ 
ing equal to the number of magnetic barriers. Let us call 
these resonances for under-barrier tunneling ‘bound-state 
tunneling resonances’ (BTRs). For the defected mag¬ 
netic barrier array, on the other hand, sharp transmission 
peaks emerge in the TGs. We call these ‘defect-induced 
tunneling resonances’ (DTRs). Due to their positioning 
within the TGs, a large modulation in transmission prob¬ 
ability is expected near their peaks. 

We next consider the effect of doping induced by ad¬ 
sorbates on the graphene sheet. In the periodic array 
case we assume that all the inter-barrier regions are ex¬ 
posed to adsorbates, while in the defected array case we 


FIG. 3. Effects of (a) the inter-barrier doping and (b) the 
local doping in the defect region on transmission probability. 
2D maps of transmission are depicted as functions of Fermi 
energy and doping potential in the case of normal incidence. 
The right insets display the doping-dependent shift of the 
transmission peaks at given Fermi energies, corresponding to 
the colored vertical lines. 

The effects of doping on the transmission probability 
are shown in Fig. [S] In the case of inter-barrier dop¬ 
ing of a periodic array without structural defects, the 
transmission spectra are entirely shifted by the doping 
potential. This shift is enough to produce strong qualita¬ 
tive changes in the transmission spectra above the barrier 
height, e.g. at the specific energy e ^ 46 meV, where the 
TG disappears and the transmission of Dirac fermions is 
almost unaffected by doping. Local doping in the defect 
region, on the other hand, does not lead to a global shift 
of the transmission spectra, but produces a clear shift of 
the DTRs: the transmission peaks are shifted, and new 
ones periodically appear as the local doping potential in¬ 
creases. 

The periodic nature of the transmission peaks is well 
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interpreted by the quantum phase through the regions 
where doping potentials are induced. The round-trip 
phase acquired by Dirac fermions while moving through 
a distance d is given by 

At a given energy, the periodicity of the potential en¬ 
ergy specific strength AU that yields Aip = 2tt can be 
found. In the particularly simple case of normal inci¬ 
dence, Eq. (IT4)) leads to a universal AU = nhvp/d where 
d = D for the inter-barrier doping case and Wd for the lo¬ 
cal doping case. Indeed, for the inter-barrier doping, we 
find AU = 41.4 meV, as verihed in the inset of Fig. [31[a), 
while for the local doping the energy interval between the 
defect-induced transmission peaks is AU = 20.7 meV, 
also shown in inset of Fig. [H^b). The same insets show 
the universality of the doping potential periodicity, being 
the same for different energies at normal incidence, albeit 
it will vary according to the incident angle. 



FIG. 4. Local-doping dependence on transmission spec¬ 
tra around different transporting gaps as a function of Fermi 
energy, in (a) the under-barrier tunneling and (b) the over¬ 
barrier tunneling regimes. The magnitude of the transmission 
probability is depicted in logarithmic scale. 

The high sensitivity of the DTRs to small local doping 
is shown in Fig. |4l in two different energy ranges cor¬ 
responding to TGs in under- (e < Et) and over-barrier 
(e > Eb) tunneling regimes. Owing to their position¬ 
ing within a TG, DTRs are sufficiently sharp, especially 
those within the lower energy TG, so that one can expect 
that the transport properties of graphene can be greatly 
modified by local doping. Note that a doping strength of 
Ud = I meV Ep = 54 meV approximately corresponds 
to a change in carrier concentration An = 8 x 10® cm“®, 
which is much smaller than the actual carrier concentra¬ 
tion of graphene ^ 2 x 10^^ cm“® at that Fermi energy. 
(See Refi^ and reference therein) 


IV. CONDUCTANCE MODULATION BY 
DOPING AT ZERO TEMPERATURE 

We studied in the previous section the transmission 
probabilities of Dirac fermions through periodic magnetic 
barrier structures in the presence or absence of struc¬ 
tural defects, and found that slight shifts of the defect- 
induced resonant peaks can produce significant changes 
to the transmission probability. In this section, we calcu¬ 
late the zero temperature ballistic conductance through 
a graphene sheet decorated with magnetic barrier arrays 
and its dependence on doping, with particular focus on 
the DTR shifts. 



0 10 20 30 40 50 60 

Fermi Energy (meV) 


FIG. 5. Ballistic conductance as a function of Fermi en¬ 
ergy through multiple magnetic structure with and without 
a structural defect. Ly = 1 fim is taken into account for 
the calculation. The inset exhibits the transmission function 
taken into account the calculation of the conductance in the 
presence of the defect. 

In the zero-temperature limit, conductance is calcu¬ 
lated from Eq. (HU, which we plot in Fig. [5] Similarly 
to the transmission spectra, defect-induced peaks that 
are modulated by local doping appear within the trans¬ 
port dips of the conductance curve. There is, however, 
an important difference: the conductance no longer drops 
to low values, especially at higher energies. At low en¬ 
ergy, on the other hand, the conductance is dominated 
by under-barrier tunneling through the barriers, resulting 
into a complete TG within which the conductance values 
almost vanish. The complete gap in low energy is due to 
the strong backscattering of the magnetic barriers, while 
the incomplete dips in high energies originate from the 
quantum interference effects in the periodic array of mag¬ 
netic barriers. The nature of the two transport regimes 
is further revealed by the symmetry of the conductance 
resonances as shown in the inset in Fig. [5] in connection 
to Eq. (fTTl) : the isotropic backscattering from the alter- 
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nating magnetic barriers in the low-energy under-barrier 
tunneling regime reflects a complete transport gap, while 
the strong dependence on the incident angle of the high- 
energy over-barrier tunneling of Dirac fermions reflects 
the incomplete drop of conductance values. 

With the possibility of a sensing application, it is nec¬ 
essary to have a figure of merit for the conductance sensi¬ 
tivity to doping. To this end, we introduce gauge factors 
{J-) as commonly used for characterizing sensing perfor¬ 
mance: 


AR/R ^ AG/G 

~ A^’ ~ A^’ 


(15) 


where R = 1/G is the graphene resistance and n = 
{Ef + Ud) /{TTfi^Vp) the local carrier concentration as a 
consequence of the local doping. These two different def¬ 
initions of the gauge factor can be selectively applied to 
specify the sensitivity, according to the way of measuring 
the transport properties of graphene. 


A. Effects of doping on the ballistic conductance 
without a defect 


on the inter-barrier doping are displayed in Fig. [6t^a). 
The conductance curves are entirely shifted as the doping 
potential varies. Figures EDb) and (c) present the gauge 
factors Tr and Tq, i.e. the sensitivity of the conduc¬ 
tance modulation by the inter-barrier doping, for various 
doping potentials. As expected, the J-q exhibits an op¬ 
posite dependence on the doping potential compared to 
Er, according to their definition (see Eq. (fT^ i. Since the 
inter-barrier doping effects lead to the entire shift of the 
conductance, very sensitive changes in the conductance 
value are obtained around the BTRs and the complete 
TG edges in the under-barrier tunneling regime. A large 
gauge factor ^ 10® is obtained near the TG edge Ep ~ 20 
meV, resulting into An = (u/Er) (AR/R) = 30 cm“^ 
for AR/R = 0.1, i.e. assuming a 10 % resistance varia¬ 
tion as the measurement resolution. Even higher gauge 
factor is obtained in the low energy edge at Ep ~ 5 meV. 
In other words, such a large value of the gauge factor 
means that the detection of the single free carrier injec¬ 
tion/extraction is possible within 3.3 mm^ area exposed 
to electrochemical adhesion for Ep = 20 mev. In the 
over-barrier regime, on the other hand, the sensitivity of 
the conductance modulation by the inter-barrier doping 
is expected to be much less sensitive, compared to the 
under-barrier regime. 



FIG. 6. (a) Ballistic conductance as a function of Fermi 

energy through periodic magnetic barriers, for different dop¬ 
ing potentials, (b) and (c) Corresponding gauge factors Er 
and Eg as a function of Fermi energy for different doping 
potentials. 

For a periodic array of magnetic barriers, the ballistic 
conductance calculated by Eq. m and its dependence 


B. Sensitivity of the conductance modulation with 
a defect: under-barrier versus over-barrier 

Next, we discuss the effects of local doping on the 
DTRs. In Figs. [TKa) and (c), the local-doping depen¬ 
dence of the conductance is shown in the under-barrier 
and over-barrier tunneling regimes, respectively, exhibit¬ 
ing shifts similarly to the transmission peaks at normal 
incidence. Despite the small change in the doping poten¬ 
tial {Ud = 1 meV), the separation of the shifted peaks are 
enough to produce a significant amount of conductance 
modulation, offering a scheme for sensitive adsorbate de¬ 
tection. 

In the under-barrier tunneling regime. Fig. [7Ka) ex¬ 
hibits a large contrast in the conductance values around 
DTRs. One can expect that the conductance abruptly 
decreases with small amount of doping because of the 
sharpness of the conductance peak in this regime. This 
ultra-sensitive conductance modulation allows us to ob¬ 
tain a very large value of the gauge factor, up to ~ 10® as 
displayed in Fig. [7Kb), but the actual magnitude of the 
conductance in this regime is too small. 

On the other hand, in the over-barrier tunneling 
regime. Fig. [TKc) exhibits that the conductance values are 
greater than those in the under-barrier regime and there 
are DTRs within the conductance dip. The DTRs are 
shifted by the local doping effects as well, but the gauge 
factor in the over-barrier regime is expected to be much 
smaller, compared to those in the under-barrier tunnel¬ 
ing regime, as shown in Fig. jT^d), due to the incomplete 
conductance drop, and thus small contrast, around the 
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FIG. 7. (a) Local-doping dependence of the ballistic conduc¬ 
tance as a function of Fermi energy through multiple magnetic 
barriers with the structural defect and (b) its corresponding 
gauge factor Tr, for under-barrier tunneling, (c) and (d) The 
same plots as (a) and (b) but for over-barrier tunneling. 

DTRs. Here, we note that for simplicity we show Tr 
only, since maximum values of J^g ars also found near 
the DTRs and as expected are almost the same as those 
of J-fl. 


C. Enhancement of the sensitivity by collimators 

The incomplete conductance drop in the over-barrier 
regime is attributed to the transport contributions com¬ 
ing from all incident angles of Dirac fermions. If these 
angular contributions could be suppressed, the transport 
dips will get deeper, resulting into a larger contrast near 
the DTRs. To this end, we propose the use of a Dirac 
fermion collimator to kill the non-zero incident angle con¬ 
tributions to the conductance. In particular, it has been 
known even a single electrostatic barrier can produce a 
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FIG. 8. (a) Effects of the collimators on the conductance as 

a function of Fermi energy. Dashed line implies the conduc¬ 
tance curve without the collimators, while solid lines indicate 
the conductance curves for different local doping potentials in 
the defect region, (b) Gauge factor Tr versus fermi energy, 
corresponding to the doping potentials. 


collimation effecij^i^. This is easily seen in the first of 
Eq. ([6]) , where for E = Uq any non-zero value of ky results 
into an imaginary kx , and thus into suppressed transmis¬ 
sion. The normal incidence transmission, on the other 
hand, remains always protected by Klein tunneling. In 
our model, we thus introduce an electrostatic barrier of 
potential height Uc and barrier width Wc = 800 nm at 
both sides of the magnetic barrier array. This leads to 
the suppression of all non-zero angular contributions to 
the conductance near Eq Uc, and thus to an increase 
of the conductance contrast around it. 

Figure [5] shows the effects of a collimator with po¬ 
tential height set to match the energy corresponding to 
a DTR. Indeed, the transmission dip becomes deeper 
in the presence of the collimator, increasing the con¬ 
ductance contrast, and thus Tr becomes large up to 
^ 10^. Such a gauge factor enables the detection of 
An = 2.24 X 10^ cm“^, for AR/R — 0.1. As experimen¬ 
tally demonstrated^, the exposure of a graphene sheet 
to 1 p.p.m. NO 2 gas results into An 5 x 10^° cm“^ 
with few % of resistance changes. With the same order 
of magnitude of AR/R, An = 2.24 x 10® cm^ is expected 
to be detected. Therefore, our results offer an oppor¬ 
tunity for sub-p.p.rn. level detection of gas molecules by 
using the doping-induced shift of DTRs. However, this 
is still much lower than the sensitivity offered by the pe¬ 
riodic array shown in the previous subsection, because of 
the still relatively large background conductance around 
the DTR, despite the use of the collimator. In closing 
this section, we remind that the application of the col¬ 
limator to the under-barrier regime would not be useful 
because complete TGs already exist there. 
































V. CONDUCTANCE MODULATION BY 
DOPING EFFECTS AT FINITE TEMPERATURE 


Our discussion regarding the conductance modulation 
has done in the zero-temperature limit. It is necessary, 
however, to also know up to what extend our results will 
be valid for the more practical cases of finite temperature. 
In this section, we investigate the finite-temperature ef¬ 
fects on the ballistic transport through the two types 
of multiple magnetic barrier structures and examine the 
conductance modulation by doping based on Eq. (HOD. 


A. Temperature effects on the conductance 
through periodic magnetic barriers 


We consider the periodic array of magnetic barriers 
with all inter-barrier regions exposed to electrochemical 
dopants (see Fig. [Ha)). 



FIG. 9. (a) Temperature dependence of the ballistic con¬ 

ductance throngh periodic magnetic barriers as a function of 
fermi energy, in the under-barrier tunneling regime, (b) and 
(c) Gauge factors Tr and To corresponding to the conduc¬ 
tance spectra for different temperatures, with Ud = 1 meV. 
Dashed curves represent Tr and Tg with Ud = — 1 meV. 

In the previous section, we found that the conductance 
modulation for the periodic magnetic barriers is expected 
to be significant around the edges of the complete TGs 
which correspond to the under-barrier tunneling. The 
calculated conductance modulation and the correspond¬ 
ing sensitivity are depicted in Fig. [HI showing their tem¬ 
perature dependence. At T = 0.1 K, large values of 


the gauge factor are observed near the TG edge, simi¬ 
larly with the results in the zero-temperature limit. This 
gauge factor peaks are due to the fact that the conduc¬ 
tance values abruptly increase near the TG edges with 
a small Ud- As temperature increases, the gauge factor 
peak is shifted to the center of the TG, exhibiting mi¬ 
nor reduction in the sensitivity up to 2 K. This shift of 
the gauge factor peak means that the variation of the 
conductance values near the TG edge becomes smoother 
as temperature increases. The sensitivity is dramatically 
reduced at T = 3 K, so that the sensitive detection of the 
presence of the adsorbates in-between magnetic barriers 
is valid for the low-temperature limit below 3 K. 


B. Temperature effects on defect-induced 
tunneling resonances 

As discussed, the key feature of highly sensitive con¬ 
ductance modulation in the defected array case results 
from the existence of DTRs and their shift by doping. 
In order to consider the possibility for finite-temperature 
gas sensing, it is necessary to see whether the DTRs sur¬ 
vive at finite temperatures. 



FIG. 10. Conductance spectra as a function of Fermi energy 
at different temperatures. Inset: detailed view of temperature 
dependence of a defect-induced conductance peak. 

Figure [10] displays the ballistic conductance through 
multiple magnetic barriers with a defect as a function 
of fermi energy for several temperatures. At low tem¬ 
perature T = 0.1 K, the distinct DTRs are clearly ob¬ 
served within the transporting dips, similar to the zero- 
temperature limit. As temperature increases, however, 
the DTRs, which are responsible for the high sensitiv¬ 
ity to the doping, become smoothed out. Thus, it is not 
expected to achieve high sensitivity with this system at 
finite temperatures. 

Indeed, the sensitivity of the conductance modulation 
near the DTR severely diminishes as temperature in¬ 
creases, as shown in Fig. [ffla). At low temperature 
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FIG. 11. (a) Temperature-dependence of the collimator ef¬ 

fects on the ballistic conductance as a function of fermi energy 
around the defect-induced peak in the over-barrier tunneling 
regime, (b) and (c) Gauge factors Tr and J-g correspond¬ 
ing to the conductance spectra for different temperatures, 
with f/d = 1 meV. Dashed curves represent J^r and J-q with 
Ud = -1 meV. 


T = 0.1 K, the doping dependence of the conductance 
still exhibits an apparent DTR within the conductance 
drop, similarly to the zero-temperature limit. The max¬ 
imum value of Fr ~ 500 gives An = 4.5 x 10^ cm“^, 
corresponding to a sensitivity reduction of a factor of 
2 compared to the zero-temperature case. Despite this 
reduction, it is still enough for sub-p.p.m. level gas de¬ 
tection as aforementioned. However, as temperature in¬ 
creases, the DTRs becomes completely smoothed out, 
eliminating any sensitivity, as seen in Fig. Illf bl and (c). 
Therefore, one can only expect highly sensitive detection 
of local doping at the low-temperature limit below 1 K. 
We note here, that the larger gauge factor values are 
achieved for negative doping potentials because of the 
asymmetric profile of the DTRs. 


C. Sensitive detection of temperature variation 



0 5 10 15 20 25 

Fermi Energy (meV) 


FIG. 12. Gauge factors (a) <I>i{ and (b) as a function of 
fermi energy with AT = 0.1 at different temperatures. 


temperature changes: 

AR/R AG/G 

= - - - = - - - 

AT/T’ AT/r’ 


(16) 


where AT is the change in temperature. Similar to Fr or 
^G^ ‘hi? and imply how sensitively the conductance is 
modulated as temperature changes. Figure [12] shows the 
calculated gauge factors as a function of Fermi energy in 
the under-barrier tunneling regime. It is clearly seen that 
<i )(3 is suitable for the temperature sensing. Large values 
of <i )(3 are achieved near the TG edge and the BTRs at 
low temperature T = 0.1 K. As temperature increases, 
the sharp and large gauge factor peaks are reduced and 
shifted to the center of the TG because the TG becomes 
narrower and shallower by temperature, as we aforemen¬ 
tioned. 

At specific temperatures, the detection of tempera¬ 
ture change is determined by and AG/G, as AT = 
(T/4>g) (AG/G). For example, at T = 0.1 K, a ~ 80 
can be found, which for 10 % conductance measure¬ 
ment resolution results into sub-mK level temperature 
detection. Obviously, the sensitivity of the tempera¬ 
ture detection becomes reduced at higher temperature, 
but the expected AT is still small below 1 K, assum¬ 
ing AG/G = 0.1. Therefore, the strong temperature- 
dependent behavior of the ballistic conductance, which 
reduces the potential for adsorbate detection, can be used 
for temperature sensing effects. 


We have found profound deductions of the gauge fac¬ 
tors at finite temperatures/ However, while this is a lim¬ 
itation for gas sensing, such a severe temperature depen¬ 
dence may become an asset which allows us to detect 
minute temperature changes. 

In order to characterize the temperature sensing abil¬ 
ity, we define the following gauge factors with respect to 


VI. SUMMARY 

This article assesses how the ballistic conductance 
through multiple magnetic barriers is modulated by lo¬ 
cal doping induced via electrochemical adsorbates on 
graphene surface. Both periodic and defected barrier ar¬ 
rays were studied. In the zero-temperature limit, large 
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sensitivity to electrochemical doping with single molecule 
detection was found at the edges of the transport gaps in 
the periodic case under inter-barrier doping, while sub- 
p.p.m sensitivity was found on the defect-induced peaks 
of the defected barrier array for local doping in the defect 
area only. In the latter case, the use of a Dirac fermion 
beam collimator (i.e. a suitable electric potential barrier) 
to suppress transport of Dirac fermions with non-zero in¬ 
cident angles was found to be necessary. 

The sensitivity of the conductance modulation is also 
discussed in the finite-temperature case. The defect- 
induced transport resonances (DTR) within the conduc¬ 
tance gaps are smoothed out as temperature increases, 
even around 1 K, and so that the sensitivity of the con¬ 
ductance modulation becomes much less than that in the 
zero-temperature limit. In the case of periodic barri¬ 
ers with inter-barrier doping in the under-barrier tun¬ 
neling regime, large values of the sensitivity gauge fac¬ 
tors can still be found below 3 K, however, as tempera¬ 
ture increases beyond 3 K, the gauge factors drastically 


decreases because of the thermal smoothing. Interest¬ 
ingly, this strong temperature dependence of the sensi¬ 
tivity is advantageous in terms of temperature sensing 
effects, allowing sub-K level temperature detection. This 
temperature sensing capability is found to increase as 
temperature decreases towards zero, pointing towards a 
temperature sensor that could be particularly useful in 
low-temperature experiments, outer space, etc. Further¬ 
more, the combination of the strong temperature depen¬ 
dence with the low thermal capacitance of the electron 
gas, may be an excellent platform for extremely sensitive 
calorimetric studies involving energy transfer between 
adsorbates, optical transitions, relaxation processes, etc, 
in low-temperature experiments. 
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Appendix A: Optical analogy of the transfer matrix 
and analytic formalism 

For plane-wave solutions, the propagating behavior of 
Dirac fermions can be governed by an optical analogy. 
The transfer matrix is also expressed by an alternative 
representation. 

As the starting point, let us consider a boundary a; = 0 
between two regions with different potentials. In the left 
and right sides of the boundary, wavefunctions are given 

by 

as well-interpreted in the main text. These plane-wave 
solutions must be matched by the wavefunction conti¬ 


nuity. The boundary condition leads to the following 
equation: 

/ 02 \ ^ 1 / Sll Sl2 ^ ^ 

\ ^2 y 2 s2 cos (j)2 \ S21 S22 y \ ) ' 

= !.«(“;), (A2) 

where 

Sll = + S2e-*‘^^ si2 = + S2e-*‘^^ 

S21 = -b 526*“^^, S22 = -b 526*"^^. 

(A3) 

The ‘interface’ matrix Ii <_2 implies scattering effects at 
the interface a; = 0 when Dirac fermions come from the 
region 2 to the region 1. Next, we look at the propagation 
of Dirac fermions through a region where potentials are 
homogeneously given. The ‘propagation’ matrix is given 

by 

/ „ikih Q \ 

Pm. = (^ 0 j . (A4) 

where li is the length of i-th region. 

Now, we can study transmission of Dirac fermions 
through a potential barrier by using the interface and 
the propagation matrix: 


03 

bs 


= P3,/3l3<-2P2,i2l2<-l 



= M 



(A5) 


The matrix M corresponds to the transfer matrix that 
describes the relation between incoming and outgoing 
waves. The matrix elements are given by 


mil = [cos (q^W) 

■ / /ss'— sine!)sin0 

-b^sm(ga;lF) 


\ COS (j) cost 


mi 2 = e ^ [cos (qxW) 


—% sin 


{q.W) ( 


ss' — sin (h sin 6 


\ cos cp cos ( 


m 2 i = —sin (gj-VF) 


ss sm m — sm ( 


m 22 = sin 


{q.W) ( 


cos (f) cos 9 
ss' sin A — sin 9 


\ cos A cos 9 


(A6) 


where all parameters qx, A^ si; S 2 are the same 
as those given in the main text. Also, W and D are the 
barrier and the inter-barrier width. 

We next formulate the transmission problem: 
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where r and t are the reflection and transmission coef- straightly seen det [P] = 1, and the determinant of the 

ficients. Solving this equation, we obtain the following interfaces matrices is obtained by 

relations between r and t: 


t = mil + mi2r, 0 = m2i + m22r, (A8) 

and we finally have: 


det [l 2 ^i] 


S COS(/) 

s' COS 9 ’ 


det [Ii^ 2 ] 


s' cos 6 

SCOS(j) 

(AlO) 


i = (A9) 

m 22 TO 22 

In order to get det [M], we take the determinant of in¬ 
dividual matrices as det [AB] = det [A] det [B]. It is 


In results, the determinant of the transfer matrix is unity, 
as expected for a lossless system. 

Thus, the reflection and transmission probabilities are 
calculated as follows: 


R=\r\^ 


T=\t\^ 


sin (qxW) {ss' sin (j) — sin 9)^ 

cos^ (j) cos^ 9 cos^ (QxW) -I- sin^ {QxW) {ss' — sin (f) sin fff' ’ 
cos^ (j) cos^ 9 

cos^ (j>cos^ 9cos^ {qxW) + sin^ {qxW) {ss' — sine/)sin0)^ ’ 


(All) 


It is easily found R + T = I, implying the flux conserva¬ 
tion. 

We formulate the transfer matrix which has the cor¬ 
respondence to the photonic analogy of multilayer struc¬ 
tures. We examine a periodic array of magnetic barriers. 
In this case, Bloch’s theorem leads to the following ex¬ 
pression: 


as 

h 




(A12) 


or more conveniently 


2 cos {KL) = mil + W 22 = Re {Tr [M]} 

= cos {kxD) cos {qxW) 

' ss' — sin (h sin 9 


- sin sin 


cos (j) cos 9 

(A15) 


where KL = kxD + ss'qxW. We now have the secular 
equation: 


mil — mi2 

W 21 m22 — 


= 0 , 


(A13) 


of which diagonalization gives the corresponding eigen¬ 
states: 


^±^KL 



This transcendental equation gives the band structure at 
the transport problem. 

We can also calculate the transmission probability 
through N barriers based on the following arguments: 



where 


(A16) 




mil sin(A^iCL) —sin[(A^—1)K’L] 
sin(/CL) 
m 2 i sin{NKL) 
sin(/CL) 


mi 2 sin(A^K’L) 
s\n{KL) 

m 22 sin{NKL)—sin[{N—l)KL] 
sin{K L) 


(A17) 


r 


In the case that the total flux is conserved, the transmis¬ 


sion probability is expressed as below: 

T 1 1 

T = 


T + R 1-h 
1 


i + I(m^)2iI 


1 + |to 21 


|2 sin‘^{NKL) 

I sin^ (iCL) 


(A18) 
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Appendix B: Band structures for magnetic 
superlattices 


As aforementioned, for an infinite number of magnetic 
barriers, i.e., magnetic superlattices, the band structures 
are obtained by 



FIG. 13. (a) Contour plot of \Re {Tr [M]}| < 2 as a function 

of energy and ky. (b)-(d) Energy dispersion versus the Bloch 
wavevector for ky = 0, 2, and 4. 


2cos{KL) = Re{Tr[M]} , (Bl) 


where K is the Bloch wavevector and L is the period of 
magnetic super lattices. The band structures of magnetic 
superlattices reflects transport properties, especially the 
nature of the transporting gap. 


Let us consider a magnetic superlattice characterized 
by the barrier width W = 50 nm, the inter-barrier dis¬ 
tance D = 50 nm, and the barrier height hypeAy/ls = 
19.8 meV. In the magnetic superlattice, the alternat¬ 
ing profile of magnetic barriers is taken into account as 
shown in Fig. [Hb). In this case, the period is given by 
L = 2(W + D). Figure [13] represents the contour plot 
of i?e{Tr[M]} as a function of Fermi energy and ky, 
and exhibits band structures as functions of the Bloch 
wavevector K. It is easily found that the band gap ~ 
2 X hypeAy/ls = 39.6 meV is consistent with the trans¬ 
porting gap shown in Fig. jSj Interestingly, there exist 
other forbidden gaps in the energy ranges [20:40] or [-40:- 
20]. Compared to the band gap between the conduction 
and valence bands, these gaps depend on ky values, (see 
Fig. [Qf bl-ldll Because of the fc^-dependence, the con¬ 
ductance values are incompletely suppressed down, form¬ 
ing transmission dips around the energy ranges [20:40] as 
depicted in Fig. [5] 

This formalism and the band structures are, of course, 
valid for only an infinite array of magnetic barriers. How¬ 
ever, if the number of magnetic barriers are large enough, 
the qualitative analysis is expected to be quite consistent 
with the finite number of barriers. 














